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Terahertz pulsed photogenerated current in microdiodes at room temperature 
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Space-charge modulation of the current in a vacuum diode under photoemission leads to the 
formation of beamlets with time periodicity corresponding to THz frequencies. We investigate the 
effect of the emitter temperature and internal space-charge forces on the formation and persistence 
of the beamlets. We find that temperature effects are most important for beam degradation at low 
values of the applied electric field, whereas at higher fields intra-beamlet space-charge forces are 
dominant. The current modulation is most robust when there is only one beamlet present in the 
diode gap at a time, corresponding to a macroscopic version of the Coulomb blockade. It is shown 
that a vacuum microdiode can operate quite well as a tunable THz oscillator at room temperature 
with an applied electric field above 10 MV/m and a diode efan of the order of 100 nanometers. 


Introduction. Space-charge limited flow in a diode has 
been a subject of investigation for over a century M. 
Nonetheless, this field of research remains quite fertile in 
terms of interesting results. The past fifteen years have 
seen considerable work on extending the classical Child- 
Langmuir law Q, e.g. for very small structures 
finite emitter area and pulse length lilI3> space-charge 
limited field emission [llMlSj. as well as development of 
novel scaling laws and investigations into the nature of 
space-charge limited flow |l3-[l7j. 

Recent research has been on the formation of regular 
electron bunches under space-charge limited conditions 
in microdiodes M- Simulations and analysis indicate 
that in diodes with gap spacing from hundreds of nm to 
a couple of microns, and applied potential of the order of 
Volts, space-charge effects can lead to the formation of 
well defined electron bunches at the cathode, that lead to 
a regularly varying anode current. Fig. [H The frequency 
of this current is in the THz band and can be tuned by 
simply varying the applied electric field inside the diode 
M- Furthermore, it appears that higher power can be 
extracted from the diode if the cathode is covered with an 
array of nanoscale emitting dots each emitting a stream 
of beamlets that is synchronized with the beamlets from 
neighboring emitters [2^ . A recent paper that presents a 
study on the maximal charge allowed for stable injection 
of electron bunches into a diode with a regular interval 
[HI has some interesting parallels to our previous work, 
although the geometry and operating parameters are not 
directly comparable. 

Although this is an interesting idea for THz generation 
and applications , there is still a fundamental issue 

that needs to be addressed. Under what circumstances 
can one expect to still see the time dependent structure 
intact? Two particular candidate causes for degradation 
of the beamlets are: velocity spread (or emittance); and 
Coulomb forces causing beamlets to expand and merge. 
In this paper we will investigate the issue by looking in 
detail at how diode gap spacing, applied potential and 
velocity spread of electrons emitted from the cathode in¬ 
fluence pulse formation and degradation. 

Physical and simulation setup. The system under 
study is an infinite parallel plate vacuum diode, with gap 



FIG. 1. The device is a vacuum diode with a disk-shaped 
emitter on an infinite cathode under a strong laser pulse. The 
electrons are extracted in bunches and are driven to the ab¬ 
sorber (anode) by the electric field E. 

spacing D, sketched in Fig.[Tl Emission from the cathode 
is restricted to a circular area of radius R = 150 nm. A 
voltage is applied across the diode gap, Vg, and it is kept 
constant throughout one simulation. The vacuum field 
for a diode without any electrons in the gap is E. 

Emission is not source limited, so we assume unlimited 
supply of electrons. This creates in turn, space-charge 
limited current in the diode gap. The flow of current is 
as follows: Once the device is turned on, a bunch starts 
to form in front of the emitter. To extract an electron 
from a particular emitter location, the surface field must 
point towards the emitter, i.e. that spot must be favor¬ 
able for emission. If so, emission is permitted and the 
emitted electron can probably be accelerated towards the 
absorber. We say probably because, the electrons already 
in the gap may position themselves in such a way that, 
later on, they can push this electron back into the emit¬ 
ter. After this, a new, random spot is checked on the 
cathode. If the new spot is favorable, as before, an elec¬ 
tron is placed 1 nm above the surface of the cathode. 
This is done to ensure that the normal component of the 
electric field on the cathode due to the newly placed elec¬ 
tron is nonzero. In reality such emitters will have surface 
irregularities which will be larger than 1 nm which per¬ 
mits us to choose such distance for electron placement. 
The numerical results are not sensitive to the exact height 
of placement as long as it is much smaller than the gap 
spacing [T9[. If an unfavorable spot is found, a new one 
is randomly picked and tested. If in hundred consecutive 
checks, no favorable spot is found, the cathode is deemed 
unfavorable for emission and the code continues to the 
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next step which is the advancement of the electrons. 

Now, the forces between all electrons in the diode gap 
are calculated together with the forces due to the applied 
vacuum field. These forces, coupled with the basic Verlet 
method [2614^ , advance the electrons for the duration of 
one time step, ts = 2.5xl0“^^s. After that, it is checked 
whether electrons have entered the anode, or if due to 
Coulomb forces they are pushed back into the cathode 
by the electrons ahead. If they did, the time and place 
of absorption is noted and the electrons are taken out of 
the system. In this moment one time step is finished, 
and a new time step begins with the emission process all 
over again. The complete explanation and setup of the 
system is given in detail in previous work, but only for 
initial velocities of electrons equal to zero [1^, [l^ . In this 
present paper we consider non-zero initial velocities. 

For the photoemission process we will divide 

the emission region into two separate regions: metal and 
vacuum. The energy 5 of electrons inside the metal cor¬ 
responds to the free electron gas and it is governed by the 
Fermi-Dirac statistics. The energy of the vacuum state 
is E'f + where Ep is the Fermi energy and 


0eff — 0w 0S 



( 1 ) 


(/)w is the work function, (ps is the Schottky work function, 
and is the height of the photoemission potential bar¬ 
rier. The energy after emission is equal to ehuj — . If 

the laser used for photoemission can be tuned such that 
huj = 5 then the electron energy will be distributed by 

the Maxwell-Boltzmann statistics jsil-lsHj. The fact that 
the Fermi-Dirac distribution is virtually identical to the 
Maxwell-Boltzmann distribution for energies higher than 
1.005F^f at room temperature justifies using Maxwell- 
Boltzmann for the electron density of states Q . A mis¬ 
match between the laser energy and the work function is 
allowable as long as it is much smaller that the thermal 
energy ksT^ ks being the Boltzmann constant and T the 
temperature. Also, the electrons available for emission 
come only from the highest energy states. The aim of 
this paper is to try to model a room temperature device, 
therefore the interest is in the high temperature tail of 
the distributions. We consider the Maxwell-Boltzmann 
distribution of electron velocities Vi at the emitter sur¬ 
face, in the spatial direction i = {x, y, z}, 


fv(Vi) 
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m being the mass of the electron. The velocities in the 
propagation direction 2 :, i. e. from left to right in Fig. 
[TJ can only have positive values, whereas those in the 
perpendicular directions x and y can take any value. 

Results and analysis. We will now show the results 
in parameter space. The temperatures tested are T = 
0,1, 2,4, 8,16, 32, 64,128, 256, 300 K. They are shown on 
the horizontal axes in the following figures. The gap spac¬ 
ing, i.e. the distance between emitter and absorber, has 


values D = 50,100, 200,400 nm, shown on the vertical 
axes. Each figure is made with a fixed value of the elec¬ 
tric field, E; = 1, 2,4,8,16,32,64,128 MV/m. 

For each combination of parameters we run six sim¬ 
ulations. The total simulated time for each case was 
^tot = 3.75 X 10“^^ s. For the systems with lowest fre¬ 
quencies (/min ~ 0.8 THz) we get ^ 300 bunching events 
and for the highest frequencies (/max ~ 2.2 THz) about 
800 bunching events in total. 
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FIG. 2. The averaged signal-to-noise ratio < SNR >. 


Averaged signal-to-noise ratio (SNR). The averaged 
SNR is a good measure of the quality of the beam 
modulation. It is calculated using the peak ampli¬ 
tude of the Fourier spectrum of the diode current 
As and the average amplitude of the noise A^, as 
SNR[dB] = 20 log]^o(^sMn)- The computational algo¬ 
rithm uses unfiltered Fourier transform and searches for 
the highest peak in a predetermined region. This re¬ 
gion can be found through the parameterization /[Hz] = 
A X (E[V/m])'^ given in [Tol where / is the frequency 
measured in Hz (with the applied vacuum field in units 
of V/m). The parameters A and a depend on the size 
of the emitter, i.e. the radius of the disk. Table [J shows 
values of these parameters for three emitter radii. If the 

TABLE I. Magnitude of the parameters A and a for different 
values of the emitter radius. 


Emitter radius [um] 

A 

a 

50 

779 X 10® 

0.539 

100 

326 X 10® 

0.580 

250 

257 X 10® 

0.575 


spectrum contains higher harmonics, as it often does, the 
algorithm chooses only the first harmonic for the analy¬ 
sis. Although higher harmonics can carry a significant 
portion of the energy, keeping with the conservative na¬ 
ture of this research, only the first, and always highest, 
peak is chosen. 

Erom Eig. [2]it can be seen that the transition from the 
parameter space where the SNR is high to that where it is 
low is more dependent on temperature (an almost vertical 
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boundary) for low values of the applied field and becomes 
more dependant on the gap spacing (an almost horizontal 
boundary) as field strength increases. The reason for this 
is that the number of electrons per bunch increases with 
applied electric field, which leads to greater intra-beam 
space-charge forces, which in turn come to dominate ther¬ 
mal effects as the cause of degradation. Although the 
definition of what constitutes a good signal-to-noise ra¬ 
tio is somewhat arbitrary and application dependant, we 
consider the SNR to be good if it is above 20 dB. 

We notice that the highest SNR for ^ = 64 MV/m 
occurs dX D = 100 nm instead of 50 nm as for other 
fields. This is because the degradation of the beamlets is 
a complex process. A long transit time allows beamlets 
to spread, but so does a greater number of electrons per 
bunch. However, since such an effect is not apparent at 
room temperature a detailed investigation is beyond the 
purpose of the present work. 
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FIG. 3. Average number of electrons per bunch < EPB >. 


Average number of eleetrons per buneh (EPB). The 
number of EPB increases with increasing electric field 
due to the fact that more electrons per bunch are needed 
to create field reversal at the emitter surface [T8|. The 
total number of electrons in the gap is proportional to 
the applied field (see Eq. (4) of [sl and Eq. (7) of [HI). 
For a constant field the total number of electrons is nearly 
constant, but a larger gap can accomodate more bunches, 
leading to fewer EPB. The number of EPB also increases 
with decreasing the temperature, since the bunches have 
more definitive structure and less smearing. Fig. [3l This 
creates robust and solid bunches which block the cathode 
for a long time before another bunch is formed. As the 
temperature is increased bunch destruction begins much 
faster after formation. Because no real blockade on the 
cathode is created another bunch forms much sooner than 
it would happen at lower temperatures. This might mean 
more bunches in the gap with increasing temperature, as 
can be seen in Fig.Hl but each bunch has fewer electrons 
and the current begins to resemble a continuous stream. 

Average number of bunehes in the gap (BIG). The 
number of BIG increases with increasing vacuum field. 
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FIG. 4. Average number of bunches in the gap < BIG >. 


gap spacing, and more subtly with temperature. As the 
vacuum field increases, the total number of electrons in 
the gap also grows and this means more bunches in the 
gap. This comes from the simple Child-Langmuir law 
0. Similarly, increasing gap spacing means more space 
and more bunches can be accommodated into it. As the 
temperature increases, the initial velocity of the electrons 
increases as well. This lowers the blocking potential from 
the charge already present in the gap, the emission is fa¬ 
cilitated, and the virtual cathode from the electrons al¬ 
ready present in the gap can be more easily overcome. If 
the temperature is low, all electrons exit the emitter with 
similar velocities and the bunch is more narrow in the di¬ 
rection of propagation. The < BIG > (average number 
per bunch) shows that the signal quality decreases with 
increasing this number. Which means that THz signal is 
the strongest at < BIG 1. 
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FIG. 5. Average peak frequency < /[Hz] >. 

Frequeney. It has been shown [l9| that the current fre¬ 
quency depends only on the electric field for zero emis¬ 
sion velocity. That is why all the plots of Fig. [5] show 
roughly homogeneous color; the frequency doesn’t change 
with changing temperature, nor gap spacing. Appar- 
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ent variations in peak frequency with increased temper¬ 
ature and gap spacing are numerical artefacts due to the 
fact that in the parameter regions where the bunches are 
poorly defined, the exact value of the peak frequency is 
hard to define. The average number of electrons in the 
diode gap (EIG) obtained for all temperatures and gap 
spacings, but at a fixed electric field, is shown in Fig. 
[6l This number increases simply due to the increasing 
space-charge limited current density given by the classi¬ 
cal Child-Langmuir law [2|. On the same figure we can 
see the average EPB. This was done as a way to cross 
check our work with previous work (sf . 

Conclusion. THz radiation with microdiodes where 
several bunches are present in the gap can be maintained 
only at temperatures below 10 K. Although such temper¬ 
atures are attainable they are impractical for real world 
operation. In this paper we showed that if the number 
of bunches in the system can be controlled to be close to 
one, the SNR significantly improves. Such systems can 


have a high SNR f> 25 dBi even at room temperature. 
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FIG. 6. Average number of electrons in the gap < EIG > (red 
circles), average number of electrons per bunch < EPB > 
(blue squares) and averaged frequency < f [THz] > (black 
triangles). Markers are slightly shifted horizontally for clarity. 
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